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На примере простой модели показано, что уравнения горения допускают решения со смещенными профилями температуры и 
участвующей в реакции горения концентрации вещества. Скорость искривленного пламени определяется однозначно, а длина 
Маркштейна оказывается существенно больше, чем в ранних теориях. Это позволяет, по меньшей мере, на порядок расширить 
область гидродинамической устойчивости по числу Рейнольдса. 


Пытаясь объяснить, почему ламинарное пламя в 
газе устойчиво к гидродинамическим возмущениям 
Дж. Маркштейн предположил, что скорость пламе- 
ни и„ зависит от кривизны К его фронта [1]. Когда 
фронт обращен выпуклой стороной к горючей сме- 
си и„ уменьшается, а если вогнутой -увеличивается. 
Это означает, что изогнутый под действием возму- 
щений фронт пламени стремится принять плоскую 
форму. В предложенной формуле Маркштейна 

и„=и,°(1 -І М К), (1) 

где и„° - скорость пламени с плоским фронтом, при- 
сутствует неопределенная величина /„ размерности 
длины (постоянная Маркштейна). Им же были сде- 
ланы первые попытки экспериментального опреде- 
ления этой величины. Дальнейшие теоретические 
поиски [2-5 и др.] приводили к выражению посто- 
янной Маркштейна в виде отношения коэффици- 
ента температуропроводности газа к на скорость и°„ 

Ім= 4 - ( 2 ) 

"и 

Но в этом случае объяснение гидродинамичес- 
кой устойчивости пламени ограничено. Экспери- 
менты показывают устойчивость пламени при чис- 
лах Рейнольдса Ке, по меньшей мере, на порядок 
превышающих критическое значение Ке сг , которое 
следует из формул (1) и (2). 

Кроме указанного Маркштейном механизма, 
приводящего к устойчивости пламени, может иметь 
место случай ухода возмущений из рассматриваемой 
области горения [3]. Дело в том, что в экспериментах 
пламя, как правило, занимает ограниченное прост- 
ранство (пламя горелки Бунзена или в трубе, свечи). 
Эти виды пламени имеют касательную составляю- 


щую скорости газа к поверхности горения, благодаря 
чему гидродинамические возмущения или затухают 
на стенках, или уходят в свободное пространство из 
области горения прежде, чем они успевают заметно 
вырасти [3]. Но такие явления наблюдается при дос- 
таточно больших числах Рейнольдса. 

Механизм ухода возмущений не объясняет пол- 
ностью наблюдающуюся гидродинамическую ус- 
тойчивость пламени. Если основываться только 
лишь на таком механизме и на формулах (1, 2), то 
мы уже при относительно небольших числах Рей- 
нольдса (Ке —10... 100) должны были видеть пламя с 
непрерывно колеблющимся фронтом. В реальнос- 
ти это не имеет место: колебание фронта, свиде- 
тельствующее о наличии гидродинамической неус- 
тойчивости, возникает внезапно и при больших 
числах Рейнольдса (Ке —ІО 3 ). 

Дальнейшие пути объяснения гидродинамичес- 
кой устойчивости пламени можно искать в самой 
теории ламинарного горения. Предложенные в ра- 
ботах [2-5] способы нахождения /„, несмотря на их 
математическую строгость расчетов, не исключают 
существования других подходов. Дело в том, что 
если в теории Зельдовича-Франк-Каменецкого од- 
номерного пламени с плоским фронтом его ско- 
рость движения определяется единственным обра- 
зом [5], то скорость распространения искривлен- 
ного пламени без привлечения дополнительного 
физического принципа не является однозначно 
определенной [6, 7]. Это обстоятельство приводит 
к неожиданному результату при исследовании не- 
одномерной диффузионно-тепловой устойчивости 
пламени [8]. Использование известной [5] схемы 
анализа на основе формулы (1) и без привлечения 
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дополнительных предположений относительно 
процессов в зоне химической реакции приводит к 
абсолютной неустойчивости пламени. 

В работе [9] было высказано предположение, что 
при искривлении фронта пламени меняется его 
структура: происходит своеобразное рассогласование 
полей температуры и концентрации (выгорания). 
Имеется в виду сдвиг распределений температуры и 
реагирующего вещества относительно друг друга и по 
направлению движения плоского исходного пламе- 
ни. Величина смещения считается пропорциональ- 
ной обратной величине энергии активации. 

Настоящая работа ставит целью дальнейшее 
развитие изложенных представлений об искрив- 
ленном пламени на основе более точного, чем в [9] 
теоретического анализа процессов во фронте неод- 
номерного пламени. 


1. Нахождение скорости пламени при слабо 
искривленном фронте. Бесконечно тонкая зона 
химической реакции 


Для описания процессов медленного горения, 
когда давление меняется слабо по сравнению с тем- 
пературой и плотностью газа, воспользуемся мо- 
делью с одностадийной химической брутто-реакци- 
ей А-^В [5]. Анализ задачи об искривленном пламе- 
ни будем проводить в линейном приближении. Это 
означает, что величина деформации фронта и вызы- 
ваемые деформацией возмущения малы, и можно 
пренебречь их квадратичными и более степенями. 


Рассмотрим модель ламинарного пламени в 
следующей формулировке, согласно которой ос- 
новное движение пламени происходит отрицатель- 
ному по направлению координаты х' [5]: 
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где Ж(А, 7) - скорость химической реакции; 
ѵ - скорость газа; Г -температура; А- концентра- 
ция реагента А; Б - коэффициент диффузии; 
0 - тепловой эффект реакции; Т - время; 
у' - координата, поперечная оси х'; с р - теплоем- 
кость при постоянном давлении. 

Примем, что химическая реакция имеет первый 
порядок: 

Ж(А,Г) = ѴѴехр(-^>> (1-2) 


где кц - константа, зависящая от свойств исходной 
горючей смеси; і:- энергия активации; Я - универ- 
сальная газовая постоянная. 


Выражения (1.1) и (1.2) в безразмерной форме 
имеют вид: 
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Символом А 0 обозначена начальная концентра- 
ция реагирующего вещества. Здесь ѵ, - масштаб 
скорости. Его детальное определение нам сейчас 
не требуется, но в дальнейших расчетах удобно 
считать ѵ, равным скорости движения стационар- 
ного пламени с плоским фронтом. Такого опреде- 
ления масштаба скорости мы будем придерживать- 
ся в дальнейшем. 

Граничные условия. Пусть пламя не ограничено 
в пространстве, а возмущения имеют локальный 
характер по направлению у. Тогда необходимые ус- 
ловия имеют следующий вид: 

х— >-оо: Ѳ = Ь = 0, х— >+сс: дѲ /дх = дЬ/дх= 0, 

у — > ісо : дѲ/ду = дЬ/ду = 0. (1-4) 


Приведем анализ искривленного фронта пла- 
мени в приближении бесконечно тонкой зоны хи- 
мической реакции и произвольных числах Льюиса. 
В этом случае достаточно рассматривать вместо 
(1.3) уравнения 


дѲ дѲ 
— — + \ѵ—~ 
ді дх 


В 2 Ѳ д 2 Ѳ ВЪ дЬ , (д 2 Ь д 2 Ь) 


(1.5) 


Наличие химической реакции учитывается раз- 
рывом (при х = 0) первых производных от темпера- 
туры Ѳ и выгорания Ъ [5]. Стационарные решения 

0° = ехр(х), Ь° = ехр(4-), х < 0; Ѳ 0 = Ъ° = 1, х > 0, ш° = 1, 
Ье 

описывают плоский фронт пламени. При малом иск- 
ривлении фронта мы ищем решение (1.5) в виде [9] 

Ѳ(Г,х,у) = ехр(х - ^) * ехр(х)(1 - §,), 

Ь(і,х,у) = сх Р(^)[і-Й- (1 ' 6) 


Здесь 4, <?2 ~ смещение профилей температуры 
и выгорания. Причем 

& = &(1+е), & = &(*, у). (1.7) 

Числовой параметр е, который считаем малым, 
но пока не конкретизируем, задает величину рассог- 
ласования распределений температуры и выгорания. 

При х— >со малые возмущения температуры и 
концентрации должны обратиться в ноль. Физи- 
чески это означает, что при малой деформации 
фронта пламени его температура в приближении 
бесконечно тонкой зоны химической реакции не 
меняется. Поэтому решения Ѳ=Ъ= 1 в области х > 0 
остаются неизменными с точностью до второго по- 
рядка по 4 и его производных. Как мы увидим ни- 
же, такая деформация приводит только к измене- 
нию нормальной скорости пламени. 
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Естественные науки 


Выражения (1.6) удовлетворяют линеаризован- 
ным граничным условиям на фронте химической 
реакции с точностью до высших порядков по Е\ и 
его производным на границе х=%,. 

Используя выражения (1.6) в ур. (1.5), учитывая 
представление \ѵ=\+ѵ/' , где \ѵ' - изменение скорос- 
ти пламени, и удерживая только линейные члены, 
приходим к следующей паре уравнений 


8^ , 8% 88, г , , 8% 

-р- + м> = — т-, + \ѵ =Ье — Нр 

81 8у 1 81 8у 2 


Принимая в расчет представления (1.7), получим 
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( 1 . 8 ) 


Если считать скорость пламени не зависящей от 
кривизны его фронта, то мы должны полагать >ѵ'=0. 
Тогда из (1.8) следует равенство <^=0, что означает 
абсолютную устойчивость пламени, т.к. его фронт 
просто недеформируемый. Полученный результат 
не зависит от того, что Ье>1 или Ье< 1 . Следова- 
тельно, <^=0 и при равенстве Ье= 1 . Значит, фронт 
пламени с бесконечно тонкой зоной химической 
реакции недеформируемый, и пламя абсолютно ус- 
тойчиво к гидродинамическим возмущениям. Этот 
вывод противоречит теории Ландау- Дарье [5, 10]. 

Пусть теперь \ѵ'^0. Преобразуя (1.8), находим 
.... (Ье-1)(1 + е) 8 2 ^ 

ѴѴ — тг. 

е ду 2 


Тогда уравнение для поверхности фронта пла- 
мени запишется как 

^О д _Х 0 , , (Ье-ЩІ + е) 

81 * 8у 2 ’ е 


Уточним вид параметра е. Его появление может 
быть обусловлено только процессами в зоне хими- 
ческой реакции, других видимых причин нет. Тогда 
при малой (порядка /К<1), но конечной толщине 
зоны химической реакции возмущения температу- 
ры и концентрации, вызванные деформацией 
фронта, будут отличаться. Величина отличия будет 
порядка р, так как это рассогласование распреде- 
лений температуры и концентрации вызваны про- 
цессами в зоне химической реакции. Тогда е~р. 
Кроме того, в предельном случае /1=0, когда е=0 и 
^=^ 2 , пламя должно быть, как мы уже видели, аб- 
солютно устойчивым. Причем, вне зависимости от 
величины отличия числа Льюиса от единицы: при 
е=0 из (1.8) следует <^,=0. Сказанное означает, что 

е = *(и-\)Р, (1.9) 

где 5 - неопределенная константа. Тогда формула для 
изменения скорости пламени приобретает форму 

ѵ , = _ 1 + Р «Ье-Щ (1Л0) 

Ре ду 2 


где индекс 1 у символа Е, опущен для простоты за- 
писи. Фронт пламени будет меняться во времени 
согласно уравнению [9] 

д 1 = п д !і пЛ ±М^ (1.11) 

81 * ду 1 ’ * Ре 


Приведенное выше предположение постоян- 
ства температуры пламени строго не обосновано. 
Это видно уже из того, что не были вовлечены в 
рассмотрение все возможные решения уравнений 
(1.5). В линейном приближении они, кроме приве- 
денных (1.6), имеют еще соответственно для перво- 
го и второго уравнений решения вида [5] 

/(х)ехр (Еіі + іку), ^-(х)ехр(О/ +іку), і = >/-Т (1.12) 

с комплексной частотой Г2 и волновым числом к, 
гдеДх), Дх) удовлетворяют обыкновенным диффе- 
ренциальным уравнениям 

^-^-(П + А: 2 )/ = 0, Ье^\-^-(^+1лк 2 )тР. 

сіх ах с/х ах 

Поэтому мы рассмотрим сейчас задачу о слабо 
искривленном фронте пламени в расширенной 
постановке. 


2. Расширенный анализ искривленного 
фронта пламени 

Пусть теперь возмущенные профили темпера- 
туры и выгорания имеют вид 

Ѳ = Ѳ\х + &) + Р(1,х,у), Ь=Ь\х+ ЕР + 0(1, х, у), 


где Ей С -новые неизвестные функции, имеющие 
смысл возмущений. Считая их малыми по абсо- 
лютной величине, представим Ѳ и Ъ в линеаризо- 
ванной форме 

Ѳ = Ѳ\х)-Е 1 сіѲ а /сІх+Г(і,х,у), 

Ь = Ь°(х)~Е 2 еіЬ 0 / сіх + С(1,х,у), 


где Ѳ°(х), Ъ\х) являются решениями уравнений 

лО ^2д0 


И , а ^- = ^\ + ]Ѵ(Ѳ\Ь°), 
сіх сіх 2 

ѵѵ° = Ре^ + \Ѵ(Ѳ° , Ь°), 


( 2 . 2 ) 


СІХ “ СІХ 2 

с первыми четырьмя граничными условиями из (1.4). 

Подставив (2.1) в уравнения (1.3) и учитывая 
равенство \ѵ=\+\ѵ' и (2.2) в промежуточных вык- 
ладках, после несложных преобразований получим 
линеаризованные уравнения 


дР , дР , сІѲ° 
д 1 дх сіх 
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где производные от скорости химической реакции 
вычисляются через стационарные распределения 
температуры Ѳ° и выгорания Ь°. 

В операторе Лапласа А. дифференцирование 
производится по переменным х, у. Выражения, 
стоящие в скобках левых частей приведенных вы- 
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ражений, следуя рассуждениям предыдущего раз- 
дела, приравняем нулю. Эта процедура даст две па- 
ры уравнений, первая пара - это уравнения для на- 
хождения смещений <^, | 2 : 


, 8% д$ 2 , 

81 ду 2 81 


■- Ье 


ду 1 


Вторая же пара - это уравнения для нахождения 
откликов температуры и выгорания на изменение 
формы пламени: 


рр+др 

81 дх 


= Д .Р + Щ-Р + 
8Ѳ° 


Ж<1 + Р* 


д\Ѵ СІЬ а е 

дь° ах ^ 


(2.3) 
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8\ѵ аѳ° 
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А, 


где учтена связь (1.7) между /(, и и опущен индекс 
"1" у смещения <*. Эволюция Е,(1, у) подчиняется 
уравнению (1.11). 


Остается теперь показать, что возмущения тем- 
пературы Р и концентрации продукта реакции С, 
во-первых, падают с течением времени при произ- 
вольных малых начальных деформациях \ плоской 
поверхности, во-вторых, они локализованы в зоне 
химической реакции, а температура пламени не 
меняется. Для этого численно решались уравнения 


8Р ! 8Р _ д 2 Р ! д 2 Р , 
81 дх дх 2 ду 2 


дК.р + Ж в + р 3 дЦ_^ (2-4) 
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Начальные условия'. 

1=0: Р=С=(), А=Ш- (2-5) 

Граничные условия для Г и С: 
х— >^х>: Р=0= 0; х->+со: дР/дх=дО/дх= 0; 
у— >+со: д Г/ ду=д С /ду= 0 . 

Граничные условия для деформации задава- 
лись в двух формах: 

1) у->+со: ді %/ду = 0; 2) у-»+со: %=0. (2.6) 

Кроме того, брались различные варианты на- 
чального отклонения 4(у). 

Стационарные распределения Ѳ°, Ь° получены 
как асимптотические (/— >аз) решения нестацио- 
нарной задачи 

^г = ^ + ЩѲ,Ь), ^- = рЛ+]Ѵ(Ѳ,Ь), (2.7) 
оі дх оі дх 


с начальными условиями 


і = 0:Ѳ=Ъ = /{х), Дх) = 


Г 1 , X — > +00 
[О, х — > -оо’ 


( 2 . 8 ) 


и граничными условиями 

х->-оо: Ѳ=Ь= 0, х->+оо: д%/дх=дЬ/дх= 0. (2.9) 

Производные от ІГв (2.4) вычисляются по ста- 
ционарным значениям температуры и выгорания и 
имеют вид 


—і 4 = 9(1-*°)— Ні ехр 

р[м+ѳ°(\-н)] 


8& 

дѴР 

8Ь° 


= -дехр 


4,11— 


Р{ н+ѳ°а-н) 


Н 1 П+Ѳ°(1 -Ц) 


Ѳ°=Ѳ°(х), Ь а = Ь\х). 


Полная схема решения поставленной задачи та- 
кова: сначала решением уравнений (2.7) с началь- 
ными (2.8) и граничными (2.9) условиями находят- 
ся установившиеся распределения Ѳ°(х), Ь\х). Пос- 
ле этого начинается решение уравнений (2.4), для 
которых начальные и граничные условия задаются 
выражениями (2.5, 2.6). 

Для численного решения всех уравнений при- 
менялась неявная схема. Причем, система (2.4) ре- 
шалась методом расщепления по пространствен- 
ной переменной Писмена-Рэдфорда [11] в квадра- 
те 40x40. Функции источников вычислялись по 
значениям Р, С на предыдущем временном слое. 

На рисунке приведены распределения Г, О в 
различные моменты времени, после возрастания в 
начальные периоды. Начальное отклонение % 0 (у) 
задавалось в виде: 

4)09=4 еслиу&[ 19,75; 20,25]; 

4)09=0, если у&[ 19,75; 20,25] ; 

Л=10 с граничными условиями 2) из (2.6). Значения Р, 
С берутся в сечении по у максимального отклонения 
А(1, у), т.е. приу=20,0, а постоянные параметры д, а, ц, 
Р полагались равными <7=20,0; 5=2,0; ^=0,15; /1=0,1. 

Как видно из приведенных графиков, возмуще- 
ния температуры Р и выгорания С сосредоточены в 
зоне пламени и убывают со временем. Изменения 
/"носят колебательный характер, причем это наб- 
людается вне зависимости от того, больше или 
меньше единицы число Льюиса. В то время как О 
(или концентрация продуктов реакции) стремится 
к нулю монотонно. Все изменения температуры и 
выгорания сосредоточены преимущественно в уз- 
кой зоне химической реакции, но незначительные 
затухания колебания могут распространяться и в 
сторону пламени, где находятся только продукты 
горения. Это означает, что природа зависимости 
скорости пламени от кривизны его фронта кроется 
в процессах, протекающих в самой зоне химичес- 
кой реакции, и не связана с распределением кон- 
центрации реагирующего вещества и температуры 
за пределами этой зоны. Поэтому утверждение 
пропорциональности с~[і вполне обоснованно. 

Обращает также на себя внимание большее по 
амплитуде значение выгорания по сравнению с 
максимумом температуры. Причем это имеет место 
вне зависимости от значения числа Льюиса. Если 
число Льюиса больше единицы, то при деформа- 
ции фронта пламени в сторону продуктов горения 
происходит уменьшение температуры и выгорания 
свежей смеси. Если же деформация происходит в 
сторону свежей смеси (рисунок), то прирост темпе- 
ратуры и выгорания положителен. 
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Рисунок. Распределения возмущений температуры Р и выго- 
рания С в различные моменты времени: а) I = 2,0 ; 
6)1=10,0. іе = 0,8 

Прежде чем делать выводы по результатам нас- 
тоящей работы, заметим, что однородные части 
уравнений для Г и С совпадают с однородными 
частями системы (2.3). Так как отклонение %({, у) 
поверхности пламени от плоского состояния быст- 
ро убывает ввиду большого значения коэффициен- 
та "диффузии" й, то условие стремления к нулю Г, 
С , в сущности, должно совпадать с условием диф- 
фузионно-тепловой устойчивости пламени, полу- 
ченным в приближении неизменности скорости 
пламени. Расчеты показали, что если числовые зна- 
чения параметров ц, /и, Р лежат за пределами облас- 
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ти устойчивого режима горения, возмущения Г я С 
с течением времени неограниченно растут. Изогну- 
тый фронт пламени при этом стремится к плоской 
форме, согласно (1.11). 

Все рассуждения проводились в рамках линей- 
ного приближения по величине кривизны фронта 
пламени, поэтому в пределах допустимости такого 
приближения решение задачи о пространственно 
искривленном ламинарном пламени может быть 
сформулировано следующим образом: 

1 . форма искривленной поверхности горения явля- 
ется решением параболического уравнения (1.11) 
с эффективным коэффициентом диффузии; 

2. скорость движения искривленного пламени яв- 
ляется суммой скорости плоского пламени и 
поправки по формуле (1.10); 

3. условие существования такого неодномерного 
пламени определяется диффузионно-тепловой 
устойчивостью ламинарного горения с посто- 
янной скоростью пламени. 

Теперь можно привести численную оценку по 
порядку величин устойчивого по отношению к 
гидродинамическим возмущениям размера Ь пла- 
мени. Для этого достаточно воспользоваться изве- 
стной формулой Маркштейна [1, 5], в которую не- 
обходимо внести множитель 1/(я/3): 

п 4к к 

' «-І5/Ч 0 ' 

Для наблюдаемого экспериментально устойчиво- 
го пламени коэффициент п теплового расширения 
находится в пределах 6. ..9, Д~0,1, к~1(Т 5 м 2 /с, 
и„°~0,1 м/с [5]. Считая п/х(п- 1)«1, находим І,~0,01 м. 
Отсюда получаем для критического числа Рейнольдса 
согласующуюся с экспериментом оценку Ке сг ~10 3 . 

Что касается определения константы 5, то этот 
вопрос остается открытым и находится на стадии 
проработки. 
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